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We study the application of the joint resummation to electroweak boson production at hadron 
colliders. The joint resummation formalism resums both threshold and transverse momentum 
corrections to the transverse momentum distribution at next-to-leading logarithmic accuracy. 
We obtain a good description of the transverse momentum distribution of Z bosons produced 
at the Tevatron collider. 

1 Introduction 

Electroweak boson production at hadron colliders serves as one of the major testing grounds 
for resummation techniques of perturbative QCD. In this talk we will describe the application 
of the recently proposed joint resummation formalism El to these type of processes and compare 
theoretical predictions of the formalism with the data on Z boson production at the Tevatron 
collider. 

It is well-known that the partonic subprocesses calculated using perturbation theory acquire 
logarithmic corrections due to gluon emission. The corrections arise from cancellations between 
virtual and real contributions at each order in perturbation theory and become large if a distri- 
bution is inspected near the phase-space boundary. The two prominent examples are threshold 
and recoil corrections, encountered in the partonic hard-scattering cross sections. The threshold 
corrections of the form a"ln^"~^(l — z)/{l — z) become large when the partonic cm. energy 
approaches the invariant mass Q of the produced boson, z = Q'^ /s — > 1. The recoil corrections, 
in turn, are of the form a" ln^"~^(Q^/(5y) and grow large if the transverse momentum carried 
by the produced boson is very small, Qt <SC Q- Thus, sufficiently close to the phase-space 
boundary, i.e. in the limit of soft and/or collinear radiation, fixed-order perturbation theory is 
bound to fail. A proper treatment of higher-order corrections in this limit requires resumma- 
tion of logarithmic corrections to all orders. Such techniques are well established both in the 



thresholdlli and in the recoil case. 

However, resummation of recoil and threshold corrections separately can lead to opposite ef- 
fects, i.e. suppression or enhancement of the partonic cross section, respectively. A full analysis 
of soft gluon effects in transverse momentum distributions dajdQ'^ dQ\ should therefore take 
these two types of corrections simultaneously into account. A joint treatment of these correc- 
tions was proposed in 011. It relies on a novel refactorization of short-distance and long-distance 
physics at fixed transverse momentum and energyGl. Similarly to standard threshold and recoil 
resummation, exponentiation of logarithmic corrections occurs in the impact parameter h space, 
Fourier-conjugated to transverse momentum Qt space, and Mellin-A^ moment space, conjugated 
to z space. This time both transforms are present, resulting in a final expression which obeys 
energy and transverse momentum conservation. Consequently, phenomenological evaluation of 
the joint resummation expressions requires providing prescriptions for inverse transforms from 
N and h spaces. This also involves specifying a border between resummed perturbation the- 
ory and the nonperturbative regime, by analyzing and parameterizing nonperturbative effects. 
Moreover, to fully define the expressions a procedure for matching between the fixed-order and 
the resummed result needs to be specified. In the following talk we will discuss these topics in 
more detail. 



The jointly resummed cross section 

In the framework of joint resummation we derive the following expression at next-to-leading 
logarithmic accuracy for electroweak annihilation §'i: 
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where cri''^ denotes the Born cross section, r = Q'^/S, and Q is the invariant mass of the produced 
boson. The flavour-diagonal Sudakov exponent E^^, at the next-to-leading logarithmic (NLL) 
accuracy in N and b was derived in§: 



E^^'iN,b,Q,t^,f^F) 



dkrp 



Aaiasikr)) In (^j + Ba{a,{kT)) 



(2) 



Dependence on the renormalization scale is implicit in Eq. (^) through the expansion of asikx) 
in powers of as{fj)- E^^ has the classic form of the Sudakov exponent in the recoil-resummed 
Qt distribution for electroweak annihilation, with the A and B functions defined as perturbative 
series in as Bi'ifl The coefficients in the expansion of these functions are the same r s in the 
pure Qt resummation and are known from comparison with fixed-order calculations t3'tiH13; at 
the NLL only the logarithmic terms with A^^\ B^^^ and A^'^^ coefficients contribute. 
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The quantity xi^^b) organizes the logarithms of and b in joint resummation i 
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x{N,b) 

where is a constant and we define 
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N = iVe'^E , 
b = bQe^^/2 



(3) 
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with "fE the Euler constant. As expected, with this choice of the form for xi^^b) the LL and 
NLL terms are correctly reproduced in the threshold limit, N ^ oo (at fixed b), and in the recoil 
limit 6 — > oo (at fixed N). 

The functions C{Q, b, N, ^, ^j-f) in Eq. (||) are chosen to correspond to the jointly resummed 
cross section for large N and arbitrary b, and to Qt resummation for b oo, N fixed: 

Ca/H{Q,b,N,fi,fiF) = ^Ca/j {N,as{n)) £jk{N,Q/x,f^F) fk/H{N,flF) ■ (7) 

They are products of parton distribution functions fk/H &t scale ^f, an evolution matrix Sji^ 
and coefficients Ca/j{N,as) with a structure of a perturbative series in Og- The first order 
expansions of the latter are given here by 



^ll. -s) = 1 + ^Cf (^-8 + vr^ + j^^^^ j = C-;- (AT, a.) , (8) 

C^^} (N,as) = —7 i r = C^}} iN.as) ■ (9) 

' 27r(A^ + l)(A + 2) 3/9 ^ ' ^ ' 

The matrix £ (iV, Q/X) A'-f) represents the evolution of the parton densities from scale ^f to 
scale Q/x up to NLL accuracy □ in In A. It is derived from NLO solutions of standard evolution 
equations 00. By incorporating full evolution of parton densities (as opposed to only the 
leading A^ part of the anomalous dimension) the cross section (|l]) correctly includes the leading 
a" ln^"'~^(A) /A collinear non-soft terms to all orders. In fact, due to our treatment of evolution, 
expansion of the resummed cross section (|l|) in the limit A — > 00, 6 = returns all 0{1/N) terms 
in agreement with the fixed-order result. Further comparison can be undertaken in the limit 
6 — > 00, A^ = when our joint resummation turns into standard Qt resummation. Consequently, 
the NLO transverse momentum distribution is recovered from the 0{as) expansion of the jointly 
resummed cross section exactly in the same way as in the Qt resummation. Outside of these 
limits, a numerical comparison between the fixed-order and the expanded jointly resummed 
expression for da/dQT at 0{as) shows, as expected, a very good agreement-especially in the 
small Qt region. 

Inverse transforms and matching 

The jointly resummed cross section (^ requires defining inverse Mellin and Fourier transforms 
so that singularities associated with the Landau pole are avoided. Similarly to threshold resum- 
mation, there are also singularities in A space associated with parton distribution functions. A 
contour for the Mellin integral in ([l|) is chosen to be analogous to the 'minimal prescription' 
contour in threshold resummation 113 

N = C + ze^"f' . (10) 

For (j) > it/2, this results in an exponentially convergent integral over A in the inverse transform, 
Eq. (|l|) for all r < 1 t3ll3. The constant C has to be such that the contour lies to the right of 
the rightmost singularity of the parton distribution functions but left of the Landau pole. 

The inverse Fourier integral from b space also suffers from the Landau singularity. We define 
this integral using the identity 

^ /'OO roo 

dHe'^^-^f{b) = 27r / dbb Jo{bQT)f{b) = vr / dbb [hi{bQT,v) + h2(,bQT,v)] f{b) , 
Jo Jo 

(11) 

and employing Cauchy's theorem to deform the integration contour along the real b axis into a 
contour in complex b plane&cl. The auxiliary functions ^1,2 (-Z; v) are related to Hankel functions 



and are defined by integrals in the complex 0-plane 
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h2{z,v) = -- I ' (^0e-*"^'°^ (12) 



The h functions distinguish between the positive and negative phases in Eq. (111]). The b integral 
can thus be written as a sum of two contour integrals (plus contributions vanishing at large \b\, 
coming from closing the contours), one integral of the corresponding integrand with hi along a 
contour C\ in the upper half of the b plane, the other integral of the integrand with /i2 along a 
contour C2 in the lower half. The Landau pole can be avoided if one defines the contours in the 
following way: 

r . h=J (0 < t < fee) r ■ h=l ^ (0 < t < 6c) 

^ • ^ ^ te-''f'>' (0 < t < cx)) ^ ■ \ be- te'^' (0 < t < 00) ' 

(13) 

where parameters be and cpi, are arbitrary as long as the contour does not run into the Landau 
pole singularity or singularities associated with the particular form (^) of the function x- 

In the joint resummation we adopt the following matching prescription between the re- 
summed and the fixed-order result: 

da da'^^ dfj^'^PC^) (^(jfixodCk) 

+ .^2.^2 ' (14) 



dQ^dQl dQ^dQl dQ^dQ^ dQ^dQ^ ' 

where da^^^ / dQ"^ dQ"^ is given in Eq. (|l]) and da'^^'^^^^ /dQ^dQj, denotes the terms resulting from 
the expansion of the resummed expression in powers of ots{^) up to the order k at which the fixed- 
order cross section do"^™'^*-'^^ ldQ'^dQ\ is taken. The above matching prescription in (iV, 6) space 
guarantees that no double counting of singular contributions occurs in the matched distribution. 

Transverse momentum distribution for Z production 

Joint resummation predictions for Z boson production compared with the latest CDF data from 
the Tevatron colliderEJ are shown in Fig. ||. Due to the contour integral prescription for perform- 
ing inverse transforms, in the framework of joint resummation one does not require any extra 
nonperturbative information to obtain predictions. This is not the case in the standard Qt re- 
summation formalism, where nonperturbative parameters are introduced to make the theoretical 
expression well defined. As shown by the dashed line in Fig. |l| the joint resummation without 
any extra nonperturbative input already provides a good description of the data, except for the 
region of very small Qtj where the nonperturbative effects are expected to play a significant role. 
However, the form of the nonperturbative input can be predicted within the joint resummation 
by taking the limit of small transverse momentum of soft radiation in the exponent, Eq. (|2|). 
Assuming moderate threshold effects the procedure gives a simple Gaussian parametrization 
FNp{b) = exp(— (76^). The value of the parameter g = O.SGeV^ is determined biL^fitting the 
predicted distribution to the data. It is very similar to the value obtained in Ref. ta, where an 
extrapolation of the Qy-resummed cross section to large b was made. The solid line in Fig. |l| 
represents predictions including the nonperturbative parametrization. In the large Qt region, 
see Fig. ^3, the joint resummation formalism with the matching prescription (|l^ ) also returns 
a very good description of data without requiring an additional switching to pure fixed-order 
result, unlike in the standard Qt resummation formalism. 

In summary, we obtain a well-defined jointly resummed cross section, valid for all nonzero 
Qt, which successfully describes the Qt distribution of Z bosons produced at the Tevatron. 
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Figure 1: CDF data0 on Z production compared to joint resummation predictions (matched to the 0{aa) 
result according to Eq. (^) without nonperturbative smearing (dashed) and with Gaussian smearing using the 
nonperturbative parameter g — 0.8 GeV'^ (sohd).The dotted line shows the fixed-order result. The normalizations 
of the curves (factor of 1.035) have been adjusted in order to give an optimal description. We use CTEQ5M 
parton distribution functions, p — pj? — Q and (j> = (j>b = 25/327r, C = 1.3, foe = 0.2/(5. 
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